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Abstract 



The neutron direct radiative capture (DRC) process is investigated, high- 
lighting the role of incident p-wave neutrons. A set of calculations is shown 
for the ^■^C(n,7) process at incoming neutron energies up to 500 keV, a cru- 
cial region for astrophysics. The cross section for neutron capture leading 
to loosely bound s, p and d orbits of ^^C is well reproduced by the DRC 
model demonstrating the feasibility of using this reaction channel to study 
the properties of nuclear wave functions on and outside the nuclear surface. 
A sensitivity analysis of the results on the neutron-nucleus interaction is per- 
formed for incident s- as well as p-waves. It turned out that the DRC cross 
section for p-wave neutrons is insensitive to this interaction, contrary to the 

case of incident s-wave neutrons. 

PACS number(s): 25.40Lw,21.10Gv,23.40.Hc 
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The direct radiative capture (DRC) process of neutrons in the keV energy region has some 
pecuharity recently revived by theoretical analysis as well as by new experimental results 
0-^. Because of the non-resonant nature of the DRC process, the complications related to 
the calculation of the compound nucleus wave function in the entrance channel are removed. 
This is a general feature of all the direct capture processes, including those induced by 
charged particle reactions. However, because of the lack of the Coulomb interaction, the 
(n, 7) reaction has salient features which makes it a unique probe for investigating nuclear 
structure information. In fact, the neutron capture process can be explored in the very 
low neutron energy region where the reaction mechanism may be fully decoupled from the 
resonance process. In this way, precise information can be obtained for the structure of the 
capturing orbit and the relative contribution of the various /-wave components to the cross 
section can be examined separately. 

Because the DRC process is essentially taking place on the nuclear surface and in the 
external region, it has been recently proposed |1|] to use this reaction channel to study the 
properties of nuclear wave functions, in connection with the discovery of the neutron halo of 
light drip-line nuclei |^ . The same kind of information can also be derived from the inverse 
reaction channel (Coulomb dissociation) where a strong enhancement of the low-lying dipole 
mode has been observed and treated as an inverse DRC process [0]. 

While the DRC of protons and alpha particles have been widely investigated in the energy 
range from a few hundreds of keV up to several MeV |^, the DRC process of neutrons has 
been mainly examined at thermal {En = 0.0253 eV) energies where s-wave neutrons are 
captured into bound p orbits. The DRC formalism for thermal neutrons has been revised 
by Raman et al. [Q. They have shown in detail how the neutron-nucleus potential strongly 
affects the capture mechanism of s-wave neutrons in light nuclei, whereas no reference to high 
energy extension nor to higher partial-wave (including p- waves) contributions was given. On 
the other hand, we may expect that, as the incoming neutron energy increases the capture 
of p-wave neutrons into bound s and d orbits comes into play and, under proper conditions, 
it can be regarded as the dominant capture process. 
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The extension of the capture models required to include p-waves and higher partial 
waves into the neutron DRC process is the main task of this note. In addition, a sensitivity 
analysis of the DRC process to the neutron-nucleus potential for energies in the keV region 
will be performed. We stress that applications of DRC models have been rarely extended to 
energies higher than thermal. In nuclear astrophysics there have been several applications 
at neutron energies of interest in the r-process nucleosynthesis P| and for inhomogeneous 
big-bang theories |jlO|- For such kind of applications it is necessary to assess quantitatively 



the DRC prescriptions in a range of energies from a few up to several hundreds of keV. 

Here we will briefly revise the DRC model for s- and p-wave neutrons and apply it to 
the calculation of the ^^C(n, 7) cross sections for transitions leading to all the four bound 
states of ^'^C . In particular we will consider realistic wave functions for the initial scattering 
state and we will focus the attention on the influence of the initial /-wave character on the 
capture cross section. 



In the early works on neutron capture reactions |]rT|-[T^ it was recognized that a capture 
mechanism, in which the incoming neutron is scattered directly into a flnal bound state 
without forming a nuclear compound state, might take place for nuclei where the flnal state 
is dominated by a strong single-particle conflguration. There have been several formulations 
of the DRC mechanism differing considerably among each other in the way the incoming 
channel and the flnal state are described In general, because the direct capture 

process is alternative to the compound nucleus (CN) formation mechanism, we can separate 
the collision matrix into two components 

U,,^f = U,^f{CN) + U,^f{DRC) (1) 

where all the quantum numbers necessary to deflne the initial and flnal states have been 
lumped into the notation i and /, respectively. The reaction cross section is given by 

^-/ = |jif^-/r (2) 

where k is the wave number of the relative motion in the entrance channel. Here we will 



deal only with the DRC part of the collision matrix. The capture cross section for emission 
of electric dipole radiation (El) in the transition i ^ f is given by 



(7, 
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k',e'\Q^A\' (3) 



where kj — ej/hc is the emitted 7-ray wave number corresponding to the 7-ray energy e-y 
and e = —eZ/A is the El effective charge for neutrons. The cross section is, therefore, 
essentially determined by the matrix elements 

gS/=<*/|T^'|^^> (4) 

where T^^ — rY'^^^O, cf)) is th electric dipole operator. Here, the initial state wave-function 
is given by a unit-flux incoming wave in the entrance channel, scattered at the origin 
by the neutron-nucleus potential. The final state wave- function ^/ is given by the residual 
nucleus (bound) final state. The radial coordinate r denotes the distance of the incoming 
neutron with respect to the target nucleus. 

The entrance channel wave function can be decomposed into spherical (Z-wave) compo- 
nents 

^i^(r) = wi(r) ^^^1^ (5) 

where Wiif) depends also on the wave number k and is written, as usual, as 

wr{T) = ^V2lTTi'[/, - IJ,0,\. (6) 

Here, the common notation for the asymptotic forms of the incoming and outgoing waves, 
respectively and O^, has been adopted 

~ exp(— i/cr H — i/vr) and Oi ~ exp(-|-iA;r — -H'k). (7) 

Ui indicates the collision matrix for the scattering process in the entrance channel, v is the 
incoming neutron velocity and k the corresponding wave number. 

The matrix elements can be decomposed into the product of three factors <5ii/ =^if ' 

Af ■ Vs. 



Radial part: Indicating with ui^{r) the radial part of the final state wave function, the 
radial overlap integral is given by 

/■oo 

lijf = Ui^{r)rwiXr)dr. (8) 

If the final state wave function is dominated by a single-particle configuration with a long 
tail outside the nuclear radius (halo), there will be a strong effect on Xij^. [||]. This will be 
shown below in the case of incoming s- and p-wave neutrons. Analytical expressions for Xi-i^ 
can be derived for specific assumptions on the initial and final state wave functions. If a 
hard-sphere model for the scattering wave function and a crude square-well model for the 
bound p orbital are assumed for the initial and final state respectively one recovers the very 



well known [13| expression for the hard-sphere capture cross section 



where y = xR- Here, x is the reciprocal attenuation length of the wave function tail 
given by X = y/2fiSn/h, where /i is the reduced mass of the system and Sn the neutron 
separation energy from the residual state. More elaborate expressions corresponding to 
different assumptions on the initial and final state wave functions, can be derived for this 



overlap integral and they can be found in the literature ||13|,|8|,p!6 . 

Angular part: For a general li If transition, the angular factor which takes into account 
the magnetic degeneracy of the final state and is averaged over the initial magnetic substates 
is given by 

Kif = ^^K^/ll^ill^^)!' = 1^(^^0101^/0)' (10) 

where (....|..) is a Clebsh-Gordan coefficient. The orbital angular momentum, I, can couple 
with the intrinsic spin, s, to give the angular momentum j = 1 + s. The corresponding jj and 
j f can be combined with the target spin I to obtain the initial and final state total angular 
momenta Jj and J/. The angular-spin coefficient for the full coupling is given by 
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\{Jf,I,jf\\Y,U,I,j,)\' _ 3 
2Ji + 1 Air 



(2/, + l)(2j- + l)(2j; + l)(2J; + l)x 



(11) 



Jf 1 J/ J/ 1 h 



where {...} are Wigner-6j coefficients. 

Final state strength: The single particle strength, S, of the final-state orbit is usually derived 
from {d,p) stripping reactions (spectroscopic factor). We note here that if the radial part 
of the matrix elements, Eq. (H), is calculated with reliable wave functions (see below) the 
spectroscopic factor can be derived from a DRC analysis of the experimental cross section. 
This technique, proposed and applied in proton capture reactions 0, could be applied in 
the neutron capture channel with the advantages already noted above. 

For thermal neutrons, the main contribution to the capture process is due to incoming 
s-wave neutrons captured into p-wave orbits and emitting El radiation. At higher neutron 
energy, however, the incoming p-wave neutrons come into play as they can be captured into 
bound s and d orbits. In ^^C both these even parity orbits are available. Moreover their 
respective levels have the characteristics of being loosely bound and with large spectroscopic 
strengths: they are good example of halo structure in excited states of stable nuclei. The 
nuclear structure information on ^^C are summarized in Table |. 

Before showing the result of the full ^^C(?t,, 7) cross section calculations we will show here 
how the s- and p-wave neutron capture is affected by the neutron-nucleus potential. 

For the purpose of the present study it is sufficient to consider single particle wave 
functions for the bound-state orbits of ^^C. They have been calculated using a Wood-Saxon 
potential with a radius parameter ro = 1.236 fm, a diffuseness d = 0.62 fm and a spin-orbit 
potential strength Vg = 7 MeV. The potential well depths were adjusted so as to reproduce 
the correct binding of the four bound states in terms of the corresponding single particle 
orbits. In the case of the 2si/2 state this gives rise to Vq = 59.23 MeV. 

In the present investigation we have treated the incoming neutron channel with the 
following approximations: 
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(a) plane wave (PW) approximation 



Ui = 1 for all / 

(b) hard-sphere (HS): scattering by an infinitely deep potential well of radius R 

Ui = e-'^'^^ for / = 

= e-"^'^^ X i±|| for / = 1 

(c) a general case in which the collision matrix is calculated numerically for a given po- 
tential of Wood-Saxon form. 

In order to see the effect of these three different treatments on the radial part of the 
matrix elements we show in Figs. |l] and |^ the calculation of the integrand of Xjj in Eq. 
(§) at En = 200 keV. Because of its dominant contribution to the integral, only the real 
(imaginary) part is shown in Fig. |l] ( |^) for s- (p-) wave capture. The cross section is 
proportional to the squared area under the curves shown in the figures. 

From Fig. |l] it is evident that for an initial s-wave state the result is strongly dependent 
on the neutron-nucleus potential adopted, namely on the different collision matrix of the 
scattering channel used in the calculation. Moreover, the behavior of the wave function 
inside the nuclear radius may result in a significant cancelation of the radial matrix element. 
This result is in full agreement with the conclusions of a detailed study of the thermal 
capture in light nuclei by Lynn et al. . We believe that this sensitivity is the main source 
for discrepancy observed between the calculated and experimental thermal cross sections 
in DRC model analysis. The infiuence of compound nucleus components in the collision 
process is, naturally, another source of uncertainty. 

On the contrary, for incoming p-wave neutrons, the radial matrix elements are essentially 
insensitive to the different collision matrices considered. This is evident, from Fig. H, for 
both the p — i> s and the p ^ d transitions. In other words, in a situation where the 
capture is dominated by a DRC process of p-wave neutrons, the cross section is not sensitive 



to the neutron-nucleus potential. Hence, the capture process is essentially determined by 
the structure of the final state wave function. In particular, the component of the wave 
function outside the nuclear radius plays the principal role in the determination of the 
capture strength. This is the main result of the present investigation. 

The results of the full calculation of the capture cross section for transitions leading to 
the four bound states in ^^C are shown in Figs. ^ and |^. 

In Fig. 1^, the capture cross sections for transitions leading to the ground state (upper 
part) and to the level at = 1.26 MeV (lower part) are shown. The curves labeled by 
(a), (b) and (c) refer to three different scattering wave functions calculated according to the 
three treatments (PW, HS and diffused Wood-Saxon potential) described above. For these 
two transitions, incoming s- and d-wave neutrons are involved in composing the initial state 
wave function. The large discrepancy between the experimental values and the calculation 
is removed only if a realistic neutron-nucleus potential is employed in the calculation. In 
our case, such an agreement was obtained when the same potential used for the bound- 
state calculation was employed (curve (c) in the figure). This is an important point because 
the obtained wave functions might have had unphysical overlaps and should have been 
orthogonalized if different potentials would have been used. 

We have extended the calculations down to thermal energy {En = 0.0253 eV) using 
the same Wood-Saxon potential and the results are shown in Table |. Considering that 
no adjustment of any parameter was performed, the results of our DRC calculation can be 
considered satisfactory. 

The relative contribution of incoming d-waves can be seen in Fig. ^. There, the dotted 
lines show the contribution of the s-wave capture component only, whereas the full line 
(c) show the contribution of both s- and rf-wave components. Thought not decisively, the 
experimental results seem to follow the increasing trend due to the onset of the higher partial 
wave component. 

Finally, the cross sections for transitions leading to the Si/2 and (^5/2 orbits are shown 
in Fig. ^ These are the transitions due to the capture from p-wave neutrons. A very 
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good agreement is found for both the reaction channels, implying reliability of the DRC 
mechanism in the energy region under consideration. In this figure the conclusion drawn 
above from the radial matrix elements calculations can be verified explicitly. The results 
of the calculations obtained using the PW, HS or the full Wood-Saxon potential are barely 
distinguishable. No major influence of the neutron-nucleus potential used is revealed by the 
calculations, very well supported by the experimental results. 

In summary, our calculations have shown that, while the DRC process of s-wave neu- 
trons is strongly biased by the incoming-neutron interaction with the target, the DRC of 
p-wave neutrons is essentially insensitive to the details of this interaction. This fact can be 
used to derive important nuclear structure information on the residual nucleus like "exotic" 
components of the neutron wave function outside the nuclear radius (i.e. neutron halo) or 
single particle strength (spectroscopic factor) of the final capturing state. 

We acknowledge many fruitful discussions with Y. Nagai and C. Coceva. This work 
has been supported in part by Grant-in-Aid for Scientific Research on Priority Areas (No. 
05243102). 
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FIGURES 




r (fm) 

FIG. 1. Real part of the integrand of Eq. (P) for incoming s-wave neutrons of 200 keV and for the 
1^3/2 orbit bound by 4.96 MeV. The scattering wave function w{r) has been calculated according 
to (a) PW, (b) HS and (c) Wood-Saxon potential prescriptions. See the text for parameters and a 
detailed explanation. 
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FIG. 2. Imaginary part of the integrand of Eq. (^) for incoming p-wave neutrons of 200 keV 
and for the 25^/2 and 1^5/2 orbits bound by 1.86 and 1.09 MeV respectively. The scattering 
wave functions w(r) have been calculated according to a PW (dashed line), HS (dotted line) 
and Wood-Saxon potential prescriptions (solid line). See the text for parameters and a detailed 
explanation. 
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FIG. 3. Neutron capture cross section of ^■^C for transitions leading to the ground state of 
""^^C {upper part) and to the level at = 3.684 MeV {lower part). The experimental values are 
from the reference Q. The curves labeled with (a), (b) and (c) correspond to the three different 
assumptions concerning the collision matrix for the scattering channel. See Figs. (|l]) and text for 
explanation. The dashed-line represents exclusively s-wave contribution to the capture. 
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FIG. 4. Neutron capture cross section of ^^C for transitions leading to the orbit {upper part) 
and to the oi'bit of ^^C (lower part). The experimental values are from Ref. |^. The dashed-, 
dotted- and solid-line correspond to the three different assumptions concerning the collision matrix 
for the scattering channel. See Fig. (|2|) and text for explanations. In this case, the capture is due 
to incoming p-wave neutrons only. 
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TABLES 

TABLE L Nuclear structure information of ^'^C and thermal neutron capture cross section of 



(MeV) 


Sn (MeV) 




Jtt 


Q a 




0.0 


4.946 


1 


1/2- 


0.77 


2.10 


3.089 


1.857 





1/2+ 


0.65 




3.684 


1.262 


1 


3/2- 


0.14 


1.10 


3.854 


1.093 


2 


5/2+ 


0.58 




Total 


3.20 


Experimental 


3.53 ± 0.07 ^ 



''From Ref. fTS] 
^From Ref. p| 
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